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It is shown that a generic black hole solution of the SU (2) Einstein- Yang-Mills (EYM) equations 
develops a new type of an infinitely oscillating behavior near the singularity. Only for certain discrete 
values of the event horizon radius exceptional solutions exist possessing an inner structure of the 
Schwarzschild or Reissner-Nordstrom type. 
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Discovered soon after the regular Bartnik-McKinnon 
(BK) solutions §, EYM black holes (BH) ||-| provided 
new insights into the BH physics related to the no-hair 
and uniqueness theorems |5j. They share sphaleronic 
properties of the BK particle-like solutions |p| and also 
exhibit an unusual discreteness ('quantization' of the YM 
field on the event horizon) due to a singular non-linear 
boundary value problem in the domain between the hori- 
zon and the asymptotically flat (AF) infinity. Still, the 
existing knowledge of the EYM BH's (unlike the BK ob- 
jects) is incomplete since only external solutions have 
been constructed so far (though a qualitative discussion 
of inner solutions is available Q). Here we present brief 
results of our investigation of the interior structure of the 
EYM BH's which reveal new surprising features due to 
coupling of non-linear fields to gravity. 

Assume the static spherically symmetric magnetic 
ansatz for the YM potential 

A = (W{r) - 1) {T^de - Tg smOdLp) , 

{T^,e are spherical projections of the SU{2) generators) 
and the following parametrization of the metric 



ds^ = {A/r^)<j^dt^ - (rVA)dr2 - r^dil^ 
where dQ^ = dO^ + sin^ Odip^, and A, a depend on 



(1) 



The field equations include a coupled system for W, A 



A{W'/ry + FW = WV/r, 
{A/r)' + 2A{W' /rf = F, 



(2) 
(3) 



where V = {W'^ — 1), F — 1 - F^/r^, and a decoupled 
equation for a: 



(Incr)' = 2W'^/r. 



(4) 



These equations admit BH solutions in the domain r > 
Th for any radius of the event horizon Vh ■ The solutions 
are specified by the number n G N oi nodes of W thus 
forming a discrete set for each r/j. Although it is not 



guaranteed a priori that the chart (|I|) is extendible to 
the full region r < r^, for AF solutions we did not meet 
any singularity in the interior region unless the genuine 
one r = is reached. In terms of coordinates one 
can find three distinct local power series solutions. The 
first one is Schwarzschild-like (S), it corresponds to the 
vacuum value of the YM field |W^(0)| — 1. Using the 
mass function m(r), A = r^ — 2mr, one gets B 



m = mo(l - 46^r^ + 8 



8fe)rV(30mo) + ©(r^), 



,V4) 



0{j 



(5) 



where mo, h are (the only) free parameters. 

The second is the Reissner-Nordstrom (RN) type of 
solution which can be found assuming the leading term 
of A to be a positive constant. This requires W{Q) = 
Wo ±1,0 and gives g 



W 

A -- 



= Wo 



- WorV(2Vo) + crV(2Wot/o) + 0{r^), 
2mor + 2(c + m^W^ /V^)r^ + 0{r^), 



(6) 



what corresponds to the RN metric of the mass mo and 
the (magnetic) charge = V^, = V{Wq). The ex- 
pansion contains three free parameters Wo, mo, c. 

We have also found the third local power series solution 
assuming a negative value for A(0) (i.e. imaginary P): 



W = Wo±r- Wor^/{2Vo) + Oir^), 
A^~Vo^T8Wor/Vo + 0{r^), 
a = ai{r^ T 4WorVl^o) + 0{r^). 



(7) 



Here there is only one free parameter (Wo) for W, A. The 
corresponding space-time near the singularity is confor- 
mal to the cylinder (after a time rescaling): 



r^{dt^ - dr^ 



de^ 



■9d^^ 



(8) 



However, one may suspect that such asymptotics can 
not correspond to a generic BH. Imposing 'boundary con- 
ditions' in the singularity we obtain the same kind of the 
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singular boundary value problem as one for external solu- 
tions where a similar role is played by the AF condition. 
The latter is known to result in the 'quantization' of the 
allowed values W"{rh)- Internal 'boundary value' prob- 
lem leads to the second 'quantization', now of the event 
horizon radius rh. Therefore, EYM BH's with S and RN 
interiors may constitute only the set of zero measure in 
the whole EYM BH solution space. For a generic rh the 
interior metric is strikingly different. 

The system (^-|^) was integrated numerically in the 
region < r < using an adaptive step size Runge- 
Kutta method for various — 10~^, 10^ starting at 
the left vicinity of the event horizon r/j with one free 
parameter Wh — W{rh) satisfying inequalities \Wh\ < 1 
and 1 — Wh < rh which are the necessary conditions for 
asymptotic flatness. For given rh, the interior solutions 
meeting the expansions may exist only for some 

appropriate Wh- A numerical strategy used to find such 
W{rh) consisted in detecting the change of sign of the 
derivative W' . In the S-case we found the curve W{rh) 
which starts at —1 asrh ^ and approaches —0.1424125 
for large rh (Fig. ||) (without loss of generality we choose 
Wh < 0). Our S-curve intersects the n = 1 member of 
the family W"^{rh), corresponding to the set of external 
AF solutions. Parameters of this BH solution are shown 
in Tab. |, its global behavior is depicted in Fig. ^. 
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FIG. 1. Wirh) for the S- and RN-type interior solutions. 
Dashed lines - W"{rh) for n — 1,2 (higher-n curves lie be- 
tween the n = 2 one and the boundary rn ~ 1 — Wu^ , dotted 
line). Note that S and RN curves W(;rh) do not merge. 




FIG. 2. The n = 1 EYM BH (S-type). 



Sigma 

Delta 


\ 




1 \ w 



FIG. 3. The n = 2 EYM BH (RN-type). 




*In later publication by Breitenlohner et al. (gr-qc/9703047) 
S-solution with n = 3 was found, so our S-solution is not 
unique as was claimed in the first version of this paper. We 
also eliminate here an irrelevant for BH solutions 'microstruc- 
ture' of the RN curve near the boundary (shown previously 
in Figs. 1,4) which was perhaps a numerical artefact. 



FIG. 4. The n = 3 EYM BH (RN-type). 

Interior solutions of the RN-type, meeting the ex- 
pansions (||) in the singularity, were found for rh > 
r* = 0.990288617. The corresponding curve W{rh) (also 
shown in Fig. ^ intersects the curves W"'{rh) for all 
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n > 2. These solutions possess an inner Cauchy hori- 
zon at some r_ < ru with |W(r_)| > 1 (Figs. ||, ||). 

Solutions of the third type (|^) were studied numeri- 
cally starting from the vicinity of the origin. The unique 
solution has been found for the horizon data subject to 
conditions \Wh\ < 1 and 1 — W"^ < for the upper 
sign m (0) and VK(0) = -0.9330656, corresponding to 
Th = 1.889088 (Fig. ^). This solution, however, does not 
meet any value W^{rh) and thus does not present a BH. 
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FIG. 5. Interior solution meeting expansions (7) at r = 0. 

Hence, we have found that the EYM BH's with the 
'standard' interiors of the S and RN types exist only for 
certain discrete values of r/j. For all other (continuously 
varying) rh one observes oscillations of A in the lower 
half-plane with an infinitely growing amplitude near the 
singularity. The oscillation region starts with an expo- 
nential fall of A which typically occurs after passing a 
local maximum r™"^ (Fig. In this regime the right 
hand side of (^) becomes comparatively small with re- 
spect to other terms, and one can simplify the system 
analytically. Neglecting the right side of (||) , one obtains 
the following (approximate) first integral of the system 

m 

Z = AVF'cr/r^ = const. (9) 

In what follows we will describe in more detail the lower- 
n case (qualitatively the behavior of A in the oscillating 
regime is the same for all n). Starting from some point 
ro close to r™""^ (ro ^ J■max^^ ^j^g quantity U = W /r 
becomes approximately constant C/ = C/q » ij^-max ^ ^ 
^^max ^ ^ Iq^qj- Then, using the Eqs. (^) and 
(^, one finds the following estimate for A valid in the 
region of the exponential fall 

A = -Crcxp[-(C/or)2] , (10) 

where C is a positive constant. From this expression it 
is clear that the fall must stop at the local minimum 

r™" = {V2\Uo\)'\ 
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FIG. 6. The beginning of A-oscillations for n — 1, rn ~ 2, 
W(rh) = -0.342072. 




FIG. 7. |A[ ^ for the second oscillation cycle of n = 1 BH 
solution with Vh = 2, Wh = -0.342072. 

According to (^, |l^ , cr decreases exponentially during 
the fall of A: 

^mm _ ^max ^_^^^^max^2^ ^ C/or™°^ > 1, 

(typically r""" < ^-^axy ^j^g^. ^mm passed, U still 
remains almost unchanged, and hence the exponential 
factor in (|l^) becomes irrelevant. It follows that after 
passing the minimum, A starts to grow linearly so that 
|A/r| « 2m is constant. This regime breaks down when 
|A| becomes comparable with V^. Then a rapid growth 
of U takes place, while a is still changing slowly. At this 
stage, according to (||), the product UA remains almost 
constant. Finally A reaches the next local maximum 
(asymptotically points of local maxima approach zero), 
and then a new oscillation period starts with a grown 
up amphtude (Tab. ||, Fig. 0). It is clear from (0) that 
a monotonically decreases exhibiting rapid falls during 
exponential falls of A and keeping almost constant values 
while A is growing up. Thus, a tends to (but does not 
reach) a zero limit as r — s- 0. 
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Although the derivative W takes rather large abso- 
lute values on some intervals, the corresponding varia- 
tion of W is still small because these intervals are also 
extremely small. All the above features (small variation 
of Z and W, constancy of U while |A| is falling down) 
become more pronounced while oscillations progress im- 
plying that both Z and W have finite limits as r — + 0. 
Then neglecting the right side of (H), omitting 1 in F, 
and replacing W by its limiting value Wq, one arrives 
at the following two-dimensional autonomous dynamical 
system 

q^p, 

p:^ {ie-" ~l)p + 2e-^'' -1/2, (11) 

where A — —{Vq/2) cxp{q), and a dot stands for deriva- 
tives with respect to r = 21n(r/i/r). This system has 
one (focal) fixed point (p = 0, q = In 2) with eigenvalues 
A — (l±iv'l5)/4, its phase portrait is shown in Fig. ^to- 
gether with an invariant set p — — ~ 1/2 correspond- 
ing to the RN-type solution. The oscillating solutions 
lie above this curve. The phase motion in this region is 
unbounded, and there are no limiting circles. The limit 
q = —oo (A = 0) can not be reached, A remains negative 
valued as r — > and passes an infinite sequence of local 
maxima and minima (Fig. Isl). 
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FIG. 8. Phase portrait of (11), RN - an invariant set corre- 
sponding to the RN-type solution (dashed - zero slope lines). 

Thus, the inner behavior of the EYM black holes shows 
the following two unexpected features. The first is 'sec- 
ond quantization' of solutions with the 'S' and 'RN' type 
interiors: such solutions exist only for discrete values of 
the horizon radius. The second is an infinitely oscillating 
nature of the generic interior metric and an associated 
new type of singularity. No Cauchy horizon is formed in 
such a regime, although the blueshift at the local min- 
ima of the mass function is extremely large and tends 



to infinity as singularity is approached. Mass function 
inflates exponentially during the first stage of each os- 
cillation cycle and after some stabilization period comes 
back to nearly zero values. The upper bound of mass 
is infinitely growing at the singularity which is space- 
like conformably with common expectations. However 
the singularity is not of the mixmaster type. It is de- 
scribed effectively by the second order dynamical system 
and hence is not chaotic, although huge metric oscilla- 
tions are encountered. 
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TABLE I. S- 


and RN-type solut 


ions. 




S-type, n = 1 


RN-type, n = 2 


RN-type, n = 3 




0.613861419 


1.273791 


1.0318420 


W{rh) 


-0.8478649145 


-0.113763994 


-0.10185163 






0.02171654 


0.08948446 


W(0) 


-1 


-1.212296124 


-1.3566052 


a(0) 


0.2263801 


5.991210 • 10"^ 


1.751928 ■ 10"=* 


Mass 


0.8807931 


1.018002 


1.000277 
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TABLE II. Oscillations parameters for ru = 2, 
Wh = -0.342072 (n = 1 BH). 







r 


W{r) 


W'ir)/r 


A(r) 




4.32 


■ 10~^ 


-0.881410 


27.89 


—2.52 • 10 ^ 


j,min 


1.12 


■ 10"^ 


-0.922056 


67.03 


-1.68 ■ 10"^ 


max 
"2 


6.63 


• 10-* 


-0.926862 


7.30 ■ 10* 


-4.16 ■ 10"" 


min 
1 2 


4.73 


■ lO"'^ 


-0.930103 


1.49 • 10'' 


-8.12 • 10**^ 


'3 


6.44 


. 10-44 


-0.930120 


4.10 • 10*^ 


-1.36 • 10"'^^ 


min 
^3 


8.81 




-0.930136 


8.01 • 10*^ 




max 
' 4 


2 ■ 10" 


1.16-10^^ 


-0.930136 


0.5 ■ 102-3210" 


-2 ■ 10-2-3210" 
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